Stress Field at a Sliding Frictional Contact: 
Experiments and Calculations 

J. Scheibert * ^ , A. Prevost, G. Debregeas, E. Katzav^ , 

M. Adda-Bedia 

Laboratoire de Physique Statistique de I'ENS, UMR 8550, CNRS/ENS/Universite 
Paris 6/Universite Paris 7, 24 rue Lhomond, 75231 Paris, France 



Abstract 

A MEMS-based sensing device is used to measure the normal and tangential stress 
fields at the base of a rough elastomer film in contact with a smooth glass cylinder 
in steady sliding. This geometry allows for a direct comparison between the stress 
profiles measured along the sliding direction and the predictions of an original exact 
bidimensional model of friction. The latter assumes Amontons' friction law, which 
implies that in steady sliding the interfacial tangential stress is equal to the normal 
stress times a pressure-independent dynamic friction coefficient /x^, but makes no 
further assumption on the normal stress field. Discrepancy between the measured 
and calculated profiles is less than 14 % over the range of loads explored. Compar- 
ison with a test model, based on the classical assumption that the normal stress 
field is unchanged upon tangential loading, shows that the exact model better repro- 
duces the experimental profiles at high loads. However, significant deviations remain 
that are not accounted for by either calculations. In that regard, the relevance of 
two other assumptions made in the calculations, namely (i) the smoothness of the 
interface and (ii) the pressure-independence of /x^ is briefly discussed. 
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1 Introduction 



The sliding contact betwee n non-conformin g elastic bodies is a classical prob - 
lem in contact mechan ics (jCattaned fll938l ) ; iMindlinI (119491 ) ; iJohnsonI fjl985l ) ; 
Hills and Nowell fll994[ )l Knowledge of the surface and subsurface stress field 
in such systems is central to solid friction, seismology, biomechanics or me- 
chanic al engineering. Typical appli cations include h a rd dis k drives (e.g. iTalke 
( Il995l )). tribological coatings (e.g. iHolmberg et al.l (119981) 1. train wheels on 
rails (e.g. Guaghano and Pau J 2007)), human jo i nts (e .gr Barbour et al.l ( 1997[) ') 



and tactile perception (e.g. Howe and Cutkosky ( 1993 ): Scheibert et al. ( 20091 )). 



Theoretically, calculations of the contact stress field in the q uasi-static stead y 
sliding regim e have been performed for both h omogeneo u s (jPoritskvl (119501 ) 



Bufier ( 1959 ): Harnilton and Goodman ( 1966h 
elastic half-spaces (iKing and O'Sullivanl (Il987l) 
(|200l [ ) \. for cvlindrical (IPoritskvl (ll95of):lBufiei] ( 
( 1966 ): 



(119661 ): 



HamiltonI (lL983h ') and l ayered 



Nowell and Hills| (1988 ): Shi and Ramalingam 



■ - - 19591) : Harnilton and GoodmanI 

King and O'SullivanI ( 1987 ): Nowell and Hills ( 1988 )), circu lar ( Hamilton and Goodman 



These calculations assume a locally valid Amontons' friction law, stating that 
everywhere within the sliding contact region, the interfacial tangential stress 
Q = f^dP with p being the interfacial normal stress and fid the dynamic friction 
coefficient. Up to now, no quantitative comparison between such calculations 
and experimental stress fields has been performed. The present work first aims 
at at filli ng this lack, by t a king advantage o f a recently p roposed experimental 
method (IScheibertI (l2008h : IScheibert et all (I2008bl . l2009[ )). which allows for di- 
rect measurements of the stress field at the rigid base of a frictional elastomer 
film. 



For such a layered system, no exact stress calculation in a steady sliding con- 
tact has been provided up to now neither. All previous works indeed rely 
on the classical Goodman's assumption which states that the n ormal dis- 
placem ents at the interface due to tangential stress are negligible (IGoodman 
(119621 )). This implies in particular that the interfacial pressure field is un- 
altered when a macroscopic tangential load is applied. For a contact be- 
tween elastic half-spaces, such a normal/tangential decoupling occurs only 
if (i) both materials are identical, (ii) both are incompressi ble or (iii) one 



of both is perfectly rigid while the other is incompressible (iBufierl ( 119591 ): 
Dundurs and Bogy ( 19691 )). For layered systems, Goodman's assumption is 
never strictly true. Howev er, it is expected to be increasingly valid (i) the 
higher the Poisson's ratio ( Kuznetsov ( 19781 )). (ii) the lower the ratio of the 
contact size a over the film thickness h or (iii) the lower the friction coefficient. 
Rigorously, one has to keep in mind that Goodman's assumption does not have 
any physical ground since it does not impose the continuity of the normal 
displacements between the two solids in contact. The present work presents 
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an exact stress analysis which, for a single linear elastic incompressible layer 
(film) under plane strain conditions, goes beyond the classical description by 
relaxing Goodman's assumption. 

In section [2], we describe the experimental setup along with the calibration 
of the apparatus. In Section [3l we present both the normal and tangential 
stress profile measurements at the base of the elastomer film obtained with a 
cylinder-on-plane contact in steady sliding. In Section HI we present the exact 
model for the quasi-static steady sliding of a rigid circular frictional indentor 
against the film. In Section the results of this exact calculation are com- 
pared to that of a semi-analytical test model implemented with Goodman's 
assumption. The measurements are directly compared to both models and 
discussed. 



2 Set-up and Calibration 



Local contact stress measurements are performed with a Micro Electro Me- 
chanical Sytem (MEMS) force sensor embedded at the rigid base of an elas- 
tomer film (Fig. [1]). The MEMS' sensitive part (Fig.[Tl inset) consists of a rigid 
cylindrical post (diameter 550 /im, length 475 fim) attached to a suspended 
circular Silicon membrane (radius 1 mm, thickness 100 /im, 330 /im below the 
MEMS top surface). When a force is applied to the post, the resulting (small) 
deformations of the membrane are measured via four couples of piezo-resistive 
gauges embedded in it and forming a Wheatstone bridge (see inset of Fig. [1]). 
The MEMS thus allows to measure simultaneously the applied stress along 
three orthogonal directions, averaged over the MEMS's millimetric extension, 
in a way that will be determined through calibration. 



In the present experiments, the MEMS sensor is located at the rigid base of a 
rough, nominally fiat elastomer film of uniform thickness h = 2 mm (~ 4 times 
larger than the post's diameter) and lateral dimensions 50x50 mm. The elas- 
tomer is a cross-linked Poly(DiMethylSiloxane) (PDMS, Sylgard 184, Dow 
Corning) of Young's modulus = 2.2±0.1 MPa and Poisson's ratio i^=0.5 
flMarklfll999h l The ratio of its loss over storage moduli, measured in a parallel 
plate rheoraeter, remains lower than ~ 0.1 for frequencies smaller than 1 kHz 
( IScheibertI (120081 )). In this range the PDMS elastomer can thus be considered 
as purely elastic. The film is obtained by pouring the cross-linker/PDMS liquid 
mix directly on the sensitive part of the MEMS (cylindrical post and mem- 
brane) so that the resulting elastic film is in intimate contact with the MEMS 
sensitive part. The parallelepipedic mold used in this process is topped with a 
Poly(MethylMethAcrylate) plate roughened by abrasion with an aqueous so- 
lution of Silicon Carbide powder (mean diameter of the grains 37 /im). After 
curing at room temperature for at least 48 hours and demolding, the resulting 
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Cylindrical glass indentor 




Fig. 1. Sketch of the experimental setup. A cylindrical glass lens (radius of curvature 
129.2 mm) to which is glued a glass cover slide is driven along the x direction against 
a rough, nominally flat PDMS elastomer film (uniform thickness h = 2 mm, lateral 
dimensions 50x50 mm) at a constant prescribed normal load P and a constant 
velocity V using a linear DC servo-motor (LTA-HS, Newport). The local normal 
and tangential stress at the rigid base of the film, respectively azz and axz-, are 
measured by a MEMS force sensor, whose sensitive part is shown in the lower inset 
(left hand), along with a sketch (right hand) showing the piezo-resistive gauges 
implementation within the Silicon membrane. P and the tangential load Q applied 
on the film are measured through the extension of two orthogonal loading cantilevers 
(normal stiffness 641ib5 N.m~^, tangential stiffness 51100ib700 N.m~^) by capacitive 
position sensors (respectively MCC30 and MCC5, Fogale nanotech). 



rms surface roughness is measured with an interferential optical profilome- 
ter (M3D, Fogale Nanotech) to be 1.82±0.10 /im. This roughness is sufficient 
to avoid any measurable pull-off force against smooth glass indentors, as dis- 
cussed in Fuller and Taborl ( 1975 ). When the film is put in contact against an 
indentor, the normal and tangential loads applied, respectively P and Q are 
measured through the extension of two orthogonal loading cantilevers (normal 
stiffness 641±5N.m~^, tangential stiffness 51100±700 N.m~^) by capacitive 
position sensors (respectively MCC30 and MCC5, Fogale nanotech). 



The stress sensing device (MEMS with its PDMS film) has been calibrated 
in an earlier work (IScheibert et al. for the normal stress only. The 
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method is recalled here and extended to the tangential stress. The surface of 
the film is indented with a rigid cylindrical rod of diameter 500 /im, under a 
normal load P. With this flat punch indentor, all sensor outputs are found to 
be linear with P. By successively varying the position of this rod along the 
X direction, and assuming homogeneity of the surface properties of the film, 
the radial profiles of the normal and tangential output voltages, respectively 
Uzz{x) and Uxz{x), are constructed point by point. These profiles are then 
compared to the results of finite elements calculations (Software Castem 2007) 
for the stress and a^z at the base of a smooth axi-symmetrical elastic film 
(with the same elastic moduli and thickness as in the experiment) perfectly 
adhering to its rigid base and submitted to a prescribed normal displacement 
over a central circular area of diameter 500 fim. For frictionless conditions, 
these numerical result s could have been obtained semi-an alytically by using 



the model developed in lFretigny and ChateauminoisI (120071 ) but finite elements 



calculations have been preferred because they allowed for variable boundary 
conditions. As expected for contact regions of dimensions smaller than the 
film thickness, the stress calculated at the base of the film are found to be 
insensitive to the frictional boundary conditions. 

The vertical dimensions of the MEMS being smaller than the thickness of the 
elastomer film, one can ignore the stress field modifications induced by the 
MEMS 3D structure and consider that the base of the film is a plane. We can 
then relate the measured output voltage U to the stress field at the base of 
the film a by writing down that 

Uaz{x, y) = AazGaz ® (Taz{x, v) (1) 

where a = x oi z. A^z and A^z are conversion constants (units of mV/Pa), Gzz 
and Gxz are normalized apparatus functions and ® is a convolution product. 
Note that we use the sign convention that azz is positive for compressive 
loading. Eqs. [T] implicitly assume decoupling between the MEMS outputs. 
This has been checked to be true for the bare sensor by submitting it to either 
a uniform pressure or a pure tangential load applied directly on the Silicon 
cylindrical post. When the MEMS is embedded in the elastomer film, this 
remains true for the normal output, as checked by applying a uniform pressure 
at the surface of the film. The analogous check for the tangential output is 
not possible because any tangential stress applied on the film surface results 
in tangential stress as well as normal stress gradients at its base, which can 
not be measured separately since they induce the same deformation mode of 
the MEMS Silicon membrane. One can still use Eqs. [T] in the limit of contact 
configurations involving small pressure gradients. This is the case when one 
uses indentors with large radius of curvature such as the cylinders considered 
in the rest of this study. In this limit, the tangential output is likely to be 
insensitive to normal stress since the Silicon sensor is much stiffer than the 
elastomer. 
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In Fourier space, Eqs. ([T]) become 



AazGaz{x, y) — J' 




(2) 



where T is the bidimensional spatial Fourier Transform, its inverse, and 
fx-, fy are the spatial frequencies in the x, y directions respectively. The 
Uzz{x, y), Uxz{x, y), (Jzz{,x, y) and axz{,x, y) fields are built from the correspond- 
ing profiles along the x axis, assuming axi-symmetry, and then transformed 
using a Fast Fourier Transform (FFT) algorithm. The rapid decay of T {azz} 
and {(Jxz} with increasing spatial frequency yields a divergence of the ratio 
in Eqs. ([2]). To circumvent this difficulty, a white noise of amplitude 10 times 
weaker than the weakest relevant spectral component is added to both terms 
of the ratio before applying the FFT. The result is found to be insensitive to 
the particular amplitude of this white noise. Azz and A^z are determined so 
that the integrals of Gzz and G^z are equal to 1. The integrals of both U^z 
and (Jxz being zero, G^z is determined up to an additive constant, which was 
taken such that G^z vanishes far from the MEMS location. 
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Fig. 2. Apparatus functions AzzGzz and A^zGxz of the MEMS sensor for (a) the 
normal stress and (b) the tangential stress, respectively. Black dots are the results 
of the calibration method. Dashed lines are the approximated apparatus functions 
used for subsequent calculations. 

Both Gzz and Gxz exhibit a bell shape with a typical width of the order of 
600 /im, comparable to the lateral dimension of the sensitive part of the MEMS 
(Fig. [2]). For subsequent calculations, Gzz is approximated by a gaussian of 
standard deviation 561 /xm (Fig. [2](a)). The shape of Gxz is more complex 
and is therefore approximated by a gaussian of standard deviation 688 /im 
decorated by a rectangular foot of lateral extent 2.7 mm and amplitude 4.1 % 
of the maximum amplitude of Gxz (Fig- El^b)). We checked that a simple 
gaussian approximation of G^jzWas not sufficient to reproduce the measured 
Uxz profile when convoluted with axz- 
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Fig. 3. Validation of the calibration procedure. Measured stress profiles under nor- 
mal loading by the rigid cylinder (P = 0.69N (•), 1.72 N (■) and 2.75 N (♦) (a) 

normal stress Szzix) and (b) tangential stress Sxz{x). Comparison is made with 
Gzz (^zzix) and Gxz "X" cr^zix) for = (solid lines) and Hs = oo (dashed lines). 
The black rectangular patches on the x axis represent the contact widths (3.00, 4.50 
and 5.34 mm for P = 0.69, 1.72 and 2.75 N, respectively) obtained from the finite 
elements calculations for = 0. 

To validate this calibration procedure, the stress profiles Szzix) = Uzz{x)/Azz 
and Sxz{,x) = Uxz{x) / A^z in the x direction for cylinder-on-plane contacts un- 
der a pure normal load are measured point by point in the same way as for the 
rod indentation. Note that the sign convention for 5*^^ is the same as for a^^. 
The indentor is a glass cover slide (thickness 150 /^m, y dimension L = 8 mm) 
glued with a very thin film of cyanoacrylate onto the cylindrical part of a 
plano-convex cylindrical glass lens of radius of curvature 129.2 mm (Fig. [T]). 
The contact length in the y direction is therefore 8 mm, a dimension which is 
large enough to create locally, at its center, a y invariant stress state, but small 
enough to make the contact insensitive to flatness imperfections at the scale 
of the elastomeric fllm lateral size. Both the glass and the PDMS surfaces are 
passivated using a vapor-phas e silanization procedure whi c h red uces and ho- 
mogenizes the surface energy f Chaudhury and Whitesides f 199ll )). Each con- 
tact is formed using the following loading sequence. The indentor is pressed 
against the PDMS fllm up to the prescribed load P within 2 % relative error. 
Due to the associated tangential displacement of the extremity of the normal 
cantilever, a signiflcant tangential load Q is induced. From this position, the 
contact is renewed by manual separation which results in a much smaller but 
flnite Q. To correct for this residual load, the indentor is displaced a few mi- 
crometers tangentially down to (5 = 0. Finite elements calculations using the 
same geometrical and loading conditions are performed with both zero and 
inflnite static friction coefficients /x^ in order to provide limiting boundary 
conditions. The calculated stress profiles azz{x) and axz{x) at the base of the 
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elastic film are then convoluted by the apparatus functions Gzz and Gxz to 
allow for comparison with the corresponding experimental measurements. The 
value = 19.00 mV/Bar obtained by deconvolution allows for the pressure 
profile measurements to lie between the /^s = and /ig = 00 limiting calcu- 
lated profiles, in the whole load range further used in this work (Fig. [3](a)). 
An equally good agreement is obtained for the tangential stress profiles with 
= 7.95mV/Bar, a value 7% higher than the one determined by decon- 
volutiord] (Fig. [3](b)). We checked that Gy^ = G^z and Ay^ = A^z- These 
apparatus functions are assumed to remain valid for contacts in the steady 
sliding regimj^. 



3 Steady sliding measurements 

The steady sliding experiments are carried out as follows. Prior to sliding, 
contacts are prepared under normal load only, ranging from 0.34 N to 2.75 N, 
using the loading sequence described in Section [2l The cylindrical indentor 
is then translated tangentially over 20 mm along the positive x direction at 
constant velocity V between 0.2mm.s~^ and 2mm.s~^. Reproducibility is such 
that Q(t) differs from less than 1 % between two successive experiments (same 
P and V). The signals display a short transient followed by a steady sliding 
regime for which both Q{t) and P{t) exhibit uncorrelated fluctuations of rel- 
ative amplitude smaller than 4 %. This observation indicates that the surface 
properties can be considered as homogeneous throughout the explored area. It 
allows us to derive the stress profiles along the sliding direction directly from 
the MEMS signals through the relation Saz{x) = UaziVt) / A^z (with a = x, 
y and z). 

Figure H] shows the measured steady sliding stress field for all 3 components 
Szz, S^z and Syz at P=1.72N and \^ = 0.4mm.s~^. They have been con- 
structed from the interpolation of 19 profiles along x at different locations 
with respect to the MEMS, with 0.5 mm steps along the y axis. Each profile 
is made of 10000 data points, one every 2 /im. The line x = corresponds to 



^ This difference is very likely due to the above mentioned fact that the MEMS' 
tangential output is sensitive to pressure gradients over the size of the sensor. These 
gradients are estimated to represent less than 6 % of the tangential output for the 
rod indentation situation used to determine A^z- For the large cylinder-on-plane 
contacts under normal loading that are considered in this calibration, the normal 
stress gradients vanish with increasing normal load. They represent at most 4 % of 
the tangential output over the whole range of P used here. 

^ In steady sliding, the normal stress gradients represent a decreasing proportion 
of the tangential output with increasing normal load, less than 16 % for P = 0.34 N, 
less than 9 % for P = 0.69 N, down to less than 4 % over P = 2.40 N. 
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Fig. 4. Measured stress field for a cylinder-on-plane contact in steady sliding regime 
at V = 0.4 mm.s~^ and P = 1.72 N. (a) Normal stress Szz (b) Tangential stress along 
the direction of movement S^z (c) Tangential stress orthogonal to the direction of 
movement Syz- Lines are iso-stress curves obtained by interpolation of 19 x-profiles 
made of 10000 data points each. The shaded zone defines the region in which quasi 
2-dimensional conditions are met. The measured field is not centered on the contact 
due to limitations in the movement of the translation stage. 

the center of the cylinder-on-plane stress profile measured under normal load, 
while the axis y = corresponds to the symmetry line of the steady-state stress 
field. These fields are to a good approximation y invariant over a width of a 
few millimeters (shaded region in Fig. H]) comparable to the extension of the 
MEMS field of integration. This observation allows us to consider that the x 
profiles at y = provide an experimental realization of a 2-dimensional (i.e. y 
invariant) cylinder-on-plane friction experiment. In the following we will focus 
on these profiles and compare them with calculated stress profiles under plane 
strain conditions. For a given P, the profiles obtained with a sliding velocity 
V in the range 0.2 < V < 2.0mm.s~^ are almost undistinguishable. Thus, in 
the following, only the profiles obtained with V = l.Omm.s"^ are shown. 

movement of the indentor 
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Fig. 5. Measured stress profiles at y = for a cylinder-on-plane contact in steady 
sliding regime &tV = 1 mm.s~^ for P = 0.34, 1.03, 1.72 and 2.75 N. (a) Normal stress 
Szz (b) Tangential stress Sxz along the direction of movement. 



Figure [5] shows the measured stress profiles Szz{x) and Sxz{,x) for 4 different 
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normal loads. For both components, the profiles exhibit a similar shape with 
a maximum at the leading edge of the moving indentor whose amplitude in- 
creases with P. The tangential component is positive throughout the contact 
whereas the normal component exhibits a negative minimum at the trailing 
edge. 



4 Exact model 
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Fig. 6. Sketch of the system considered in the exact model. An elastic film is perfectly 
adhering on its rigid base {z = 0). At its surface {z = h) it is stress free outside of 
the contact region (|a;| < ah), with a being a result of the calculation. Within the 
contact region, the normal displacements Uz are prescribed and in steady sliding 
(^xz = l^d<^zz is assumed everywhere at the interface, being the dynamic friction 
coefficient. 

To allow for a direct quantitative comparison with the previous experimental 
stress profiles we have developed the following bidimensional exact model (Fig. 
E]). A linear incompressible elastic film, of thickness h and Young's modulus 
is loaded under plane strain conditions by a rigid circular body of radius 
R moving at a constant velocity V . We postulate quasi-static motion, i.e. 
the characteristic time h/c for sound waves of velocity c to travel across the 
film is assumed to be smaller than the characteristic time a/V associated 
with the indentor motion, so that the elastic film is at equilibrium at all 
times. The problem is made dimensionless by expressing the coordinates (x, z), 
displacements Ui{x,z) and stress aij{x,z) in units of h, h^/2R and Eh/6R, 
respectively. 



The constitutive equations for the elastic film can be written as 

dui duj 



dxj dxi 



(3) 
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where E is the pressure. The equihbrium equations in the film and the condi- 
tion of incompressibility are 



VS = Am (4) 

V-M = (5) 

We specify the following boundary conditions 

M,(x,0) = n,(x,0) = (6) 

Cr^^(x, 1) + fld(Tzz{x, 1) = (7) 

(Tzz{x, 1) = for |x| > a (8) 

Uz{x, 1) = —u{x) for \x\ < a (9) 



where z = and z = 1 correspond to the locations of the base and the surface 
of the elastic film respectively. Eq. accounts for the perfect adhesion of 
the film to its rigid base, Eq. ([7]) corresponds to Amontons' law of friction 
with a dynamic friction coefficient fid, Eq. ([8]) insures that the surface of the 
film is traction-free outside the contact zone and Eq. ([9]) defines the normal 
displacement induced by the indentor over the contact zone of width 2a. For 
a circular rigid indentor the normal displacement has a parabolic profile given 
by 

u{x) = - - {x - Xof (10) 
a 

where xq represents the asymmetry of the steady sliding contact and a = 
h'^/2R6 with 6 being the normal displacement of the indentor. Both xq and a 
are selected by the system for a given width of the contact zone area a and 
friction coefficient /x^. 



As suggested by the strip geometry and the boundary c onditions, the resolu- 

tion in volves the use of Fourier sine and cosine transforms (lAdda-Bedia and Ben Amar 
(I2OOII )). Any spatial distribution function D{x,y) of the problem (displace- 
ment, strain or stress) may be decomposed into 



D{x,y)= r D^^\k,y) cos kxdk+ H D^'\k,y) sin kxdk (11) 
JO Jo 

Substituting this representation into the bulk equations Eqs. ([2])^® and the 
boundary conditions ([6])-([9]) and exploiting the parity properties of the sine 
and cosine functions, lead to the following equations 

r ai'^Xk A) cos kxdk = \x\ > a (12) 
Jo 

r ai'Xk A) sin kxdk = \x\ > a (13) 
Jo 
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and 



= -\ Hx)+u{-x)] 
-f,,F,{k)a^f]{kA) + Fo{k)a^J{k,l) 
= — - [u{x) — u{—x)] 



coskx 
Ixl < a 



sin kx 



dk 



2k 

\x\ < a 



(14) 



(15) 



where 



Foik) 

m) 



sinh(2A;) - 2k 
cosh(2A:) + 1 + 2A;2 
2^ 

cosh(2A:) + 1 + 2/^2 



(16) 
(17) 



The conditions ffT^ [T^ are identically satisfied by 

a^:}{k,l)= r(P{t)Mkt)dt (18) 
Jo 

ail\k,l) = r tilj(t)Ji(kt)dt (19) 
Jo 

irrespective of (j){t) and ip{t), with Jo(3^) and Ji(x) being the Bessel functions 
of the first kind. The functions (j){t) and 'ilj{t) now become the unknowns in 
the problem. 



In two-dimensional contact problems, the indentation depth is undeterminate, 
which requires differentiating the boundary conditions ffni) - f[T5]l with respect 
to X before replacement into the representation ffT8D-(fT9l) . One t hen classically 



gets a set of coupled integral equations (see e.g. Spencd (j 19751 ) and Gladwelll 



(jl980l )). that are here of Abel type which fix the functions 0(t) and ip{t). 



Inverting this set of equations using the Abel transform yields 

0(x) + r Moo(x, t)0(t) dt + lid r Mio{x, t)ip{t)dt = -4x (20) 
Jo Jo 

ij(x)+Hd[ Moi(x,t)(j)(t)dt- [ Mii(x,t)'ip(t)dt = (21) 
JO Jo 



where Mij{x,t) = (-l)V--'f /o°° ^(^i-il(^) - Sij)Ji{kt)Jj{kx)dk. Eqs. ((20])- 
(12T!) are independent of the parameters xq and a which allows to solve them 
once the constant fid and a are fixed. This simplifies the numerical scheme. 
Then, xq and a are fixed a posteriori by using Eq. [H] and the derivative of 
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Eq. [15] with respect to x at, say x = 0. This leads to the following equations 

ra rco ra roc 

Axo=i^d Ht) Fi{k)Jo{kt)dkdt- ttpit) Fo{k)Ji{kt)dk dt (22) 
Jo Jo Jo Jo 

-=xi- r (t){t) r ^^Jo(kt)dkdt- fid rtijit) r ^^^Jiikmkdt 

a Jo Jo 2k Jo Jo 2k 

(23) 



The displacement and stress fields can be easily expressed as functions of 0(x), 
ip{x), xq and a and thus can also be calculated numerically. The lineic normal 
load Pl applied to the film surface can then be calculated using the following 
expression 

a 

Pl = - f_ <yzz{x, 1) dx = -K j (j){t)dt. (24) 

" 

Using the constitutive equations and providing simple algebraic transforma- 
tions the normal stress 0"^^(x, 0) and the tangential stress cr^zix, 0) at the rigid 
base are given by 



(^zzix, 0) = y [Zi{x, t) + fidZsix, t)] (l){t)dt - j [fidZ2ix, t) - Z^{x, t)\ tip{t)dt 



(25) 

a a 

0-xzix, 0) = - y" [f^dZ^ix, t) + Z3{x, t)] (f){t)dt + J [Z2{x, t) - fidZeix, t)] tip{t)dt 



(26) 

where the kernels Zi{x,t) are explicitly 

oo 

Zi{x,t) = J A{k)cos{kx)Jo{kt)dk, (27) 



oo 

Z2(x,t) = J B{k)cos{kx)Ji{kt)dk, (28) 



oo 

Z3{x,t) = J B{k)sm{kx)Jo{kt)dk, (29) 



oo 

Z4(x,t) = J A{k)sm{kx)Ji{kt)dk, (30) 



oo 

Z5{x,t) = J C{k)cos{kx)Jo{kt)dk, (31) 



oo 

Zeix,t) = J C{k)sm{kx)Ji{kt)dk, (32) 
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with A{k), B{k) and C{k) being 



m 



2(cosh(/c) + /csinh(A;)) 
cosh(2A;) + 1 + 2fc2 



(33) 



2/;; cosh(/c) 



(34) 



C{k) 



cosh(2A:) + 1 + 2A;2' 
2(cosh(/c) — ksmh{k)) 



(35) 



cosh(2A:) + 1 + 2A;2 



In practice, the input parameters of the model are chosen to be /i^ and P^, 
and the resulting normal and tangential stress profiles at the base of the film 
are derived. 



5 Discussion 

We recall here that the calculation presented in the previous section is the first 
one relaxing Goodman's assumption for the frictional steady sliding of a lay- 
ered material. In order to assess the impact of this increment on the mechanical 
description of such contacts, we directly compare, for various combinations of 
the input parameters /x^ and Pl, the stress profiles obtained from both our 
exact calculation and an additional calculation derived along the same lines as 
the exact one but with Goodman's assumption. The latter test model, referred 
to as Goodman's model is detailed in Appendix \^ 

Figure [7| shows the normal stress profiles crf^ = (Tzz{x, 1) at the surface of 
the film. For each normal stress profile, the corresponding tangential stress is 
obtained by multiplying the former by the friction coefficient /i^, i.e. a^.^ = 
^dO'lz following Amontons' law - see Eq. [71 As expected, for /i^ = 0, the exact 
calculation matches Goodman's result and yields symmetric fields with an 
integral (area below the curve) equal to Pl- For increasing /i^ at constant Pl, 
the profiles maintain their integral while becoming increasingly asymmetric, 
with a growing maximum shifting towards the leading edge of the moving 
indentor. A similar behavior for the envelope is observed for an increasing Pl 
at constant fid- Interestingly, Goodman's model deviates significantly from the 
exact one, even in the favourable situation considered here where the material 
is incompressible and the film is relatively thick. 

Figure [HI shows both the normal and tangential stress profiles, azz{x,0) and 
(Jxz{x, 0), at the base of the film, where the the stress a is actually measured. 
a is related to at the free surface of the film through a convolution with 
the Green function for an elastic membrane of thickness h. Since the latter 
has a typical width c::^ h, a cannot exhibit spatial modulations over length 
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Fig. 7. Normal stress profiles a^.^ = ^^^(x, 1) at the surface of the film, calcu- 
lated with the exact model (solid lines) or with Goodman's model (dashed lines), 
(a) lid increases from 0.3 to 3.0 with steps of 0.3 for the same lineic normal load 
Pl = 200 N.m~^ . For all cases, the contact radius is 2.36 it 0.03 mm. (b) Pl increases 
from 20 to 380Pa.m~^ with steps of 40Pa.m~^ for the same friction coefficient 
/id = 2.0. For the exact model, contact widths are 2.60, 3.40, 3.88, 4.24, 4.52, 4.80, 
5.02, 5.24, 5.44 and 5.62 mm respectively. For Goodman's model, contact widths 
are 2.68, 3.52, 4.00, 4.36, 4.66, 4.90, 5.12, 5.32, 5.52 and 5.68 mm respectively. For 
all these graphs, the following parameters were used: E = 2.2 MPa, R = 130 mm, 
h = 2 mm. 
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Fig. 8. (a) Normal stress profiles fJ2^(x,0) and (b) tangential stress profiles along 
the direction of movement axz{.x,Qi) calculated at the base of the elastic film with 
the exact model (solid lines) or with Goodman's model (dashed lines), increases 
from to 3.0 with steps of 0.6 for the same lineic normal load Pl = 200 N.m~^. The 
contacts widths are equal to that given in the legend of Fig. [7|^a). The following 
parameters were used: E = 2.2 MPa, R = 130 mm, h = 2 mm. 

scales smaller than h= 2 mm. The spatial resolution of the MEMS (~ 1 mm) 
is therefore sufficient to probe the stress field a at the base of the elastic film. 
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Fig. 9. (a) Normal stress profiles crzz{x,0) and (b) tangential stress profiles along 
the direction of movement axzix,0) calculated at the base of the elastic film with 
the exact model (solid lines) or with Goodman's model (dashed lines). Pl increases 
from 20 to 380Pa.m~^ with steps of 80Pa.m~^ for the same friction coefficient 
Hd = 2.0. The contact widths are equal to that given in the legend of Fig. [TJb). The 
following parameters were used: E = 2.2 MPa, R = 130 mm, h = 2 mm. 



For /irf = 0, the normal stress profile is symmetric with an integral equal to Pl 
whereas the tangential stress profile is antisymmetric with a vanishing integral. 
For a given lineic load Pl, an increasing fid qualitatively results in growing 
additional contributions to the profiles, anti-symmetric for the normal stress 
and symmetric for the tangential stress. The integral of the normal stress 
profile remains equal to Pl while the integral of the tangential stress profile 
becomes HdPi- Similar features are observed on Fig. [9], which shows azz{x, 0) 
and (Jxz(x,0) for an increasing lineic load Pl and a given friction coefficient 
fid- Goodman's model yields qualitatively similar results but with growing 
errors for increasing Pl or fid- In particular, Goodman's model underestimates 
the amplitude of the maxima of both stress components at positive x and 
overestimates the amplitude of both the negative part of the normal stress 
and the dip in the tangential profiles at negative x. 

The measured stress profiles Szz{x) and Sxz{x) along ?/ = can be now quan- 
titatively compared to the stress profiles (Jzz{x) and (7xz{x) calculated at the 
base of the elastic film and convoluted with the apparatus functions Gzz and 
Gxz determined in Section [2l In the limit of a bidimensional geometry, the 
input parameters used in the calculation - namely the applied lineic load Pl 
and the dynamic friction coefficient fid - should ideally be deduced from the 
macroscopic measurements of P (the normal load) and Q (the tangential load) 
by using P/ L and Q / P respectively, with L being the contact length. This ap- 
proach yields inconsistent stress profiles for two reasons. First, with our finite 
sized punch experimental system, the contribution of edge effects to the total 
normal load P is not negligible. For a given x, the interfacial pressure has a 
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minimum around y = 0, so that P/L over-estimates the effective hneic load 
at the location of the measured profile. Second, the measured macroscopic 
friction coefficient Q/P turns out to be a decreasing fonction of P (and thus 
of the local pressure), assuming values from 1.5±0.1 at P = 0.34N down to 
1.36±0.04 at P = 2.75N, which a r e typi c al for PDMS on glas s stead y sliding 
contacts (see e.g. Galhano et al. ( 2003 ): Wu-Bavouzet et al. ( 200?! )). These 
averaged values under-estimate the effective friction coefficient at the location 
of the measured profile since the pressure has a minimum around y = 0. To 
circumvent this difficulty, we extracted P^ and /i^ from the measured stress 
profiles as Pl = S^zdx and /i^ = S^zdx/ S^zdx. With such defini- 
tions, Pl is found to increase from 20 to 220 N.m~^ and fid to decrease from 
2.6 to 2.0 when P varies from 0.34 to 2.75 A^. 




Fig. 10. Measured stress profiles (o, for clarity only one percent of the data points 
is shown) at y = (a) Szz and (b) Sxz in steady sliding regime for increasing normal 
loads (P = 0.34N in black, 1.03 N in red, 1.72 N in blue and 2.75 N in green) and 
V = 1 mm.s~^. Comparison is made with (a) Gzz ^ ctzz and (b) Gxz ® ctxz where azz 
and Gxz are computed from the exact model (solid lines) or from Goodman's model 
(dashed lines). 

Figure [10] shows the measured profiles together with the predicted stress pro- 
files convoluted with the apparatus functions, for both our exact model and 
Goodman's model. The two calculations predict profiles in reasonable agree- 
ment with the experimental ones. In particular, they account for both the neg- 
ative part of Szz{,x) and the dip of Sxz{x) at negative x. In order to quantify 
the deviations between the experimental and calculated profiles, we compute 
the quantity x - 



T,i{Ei — Ciy/TiiEf, where Ei are the experimental data 
points and Ci are the calculated ones. For the tangential stress, both models 
yield similar values of x=ll i 1%, with no clear load dependance. For the nor- 
mal stress profiles, the exact model yields an almost constant x=ll ± 3 % over 
the range of normal loads P explored. For Goodman's model, x increases with 
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the load, between 12 and 28 %, indicating a decreasingly good fit to the ex- 
perimental data with increasing P. The exact model is therefore the one that 
follows most closely the evolution of the experimental profiles with increasing 
normal load (Fig. [TUi) . which is consistent with the fact that Goodman's as- 
sumption is expected to fail as the ratio of contact size a to film thickness h 
becomes large. 



Although the exact model accounts for the data better, non-negligible ro- 
bust deviations are observed for which we do not have any definitive expla- 
nation. Two central assumptions used in both models are however amenable 
to refinement and may explain the observed deviations. First, the interface 
is assumed to be molecularly smooth whereas the surface of the elastomer 
exhibits a micrometric roughness. The resulting multicontact interface is thus 
expected to exhibit finite compressive and shear compliances. This feature 

has been shown to modify, w i th respect to srnooth c o ntacts, both the stress 

( Greenwood and Tripp ( 1967[): Scheibert et al. J 2008b ): Chateauminois and Fretigny 
( 20081 )) and displacement ( Scheibert et al. (2008a)) fields. These effects are ex- 
pected to induce vanishing corrections at increasingly high loads. The second 
questionable assumption is the existence of a single pressure-independent fric- 
tion coefficient. This is clearly at odds with the observed decrease oi Q/P as 
a function of P. Such a be havior is usually attributed to t he finite adhesion 
energy of the interface (e.g. lCarbone and Mangialardil (120041 )). and is sensitive 
to the geometrical properties of the film roughness. 



6 Conclusion 



This work provides the first spatially resolved direct measurement of the stress 
field at a sliding contact. The choice of a cylinder-on-plane geometry has al- 
lowed us to quantitatively compare the profiles measured at the center line of 
the contact with bidimensional calculations. An exact model was developed 
to predict the stress field at the sliding contact assuming linear elasticity and 
a locally valid Amontons' friction law, but without the classical Goodman's 
assumption on the normal displacements. This model correctly captures the 
measured stress profiles with typical deviations of less than 14%. In the range 
of loads explored experimentally, this calculation does not differ drastically 
from the classical calculation involving Goodman's assumption. However, the 
present model is expected to provide significant improvements over Good- 
man's model as the thickness of the film is further reduced or as the load is 
further increased. In these cases, Goodman's assumption becomes increasingly 
inaccurate. 

Robust deviations between the experiments and the model have been briefly 
discussed along two lines, namely the finite compliance of the multicontact 
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interface and the pressure-dependence of the friction coefficient. However, the 
cyhnder-on-plane experiment described here, which was specifically designed 
to allow for a comparison with bidimensional models, is not best suited to 
study such fine effects. As discussed, the resulting edge effects do not allow 
one to use well-controlled or measured macroscopic quantities, e.g. P and Q, as 
input parameters in the models. This could be done for instance with a sphere- 
on-plane geometry, but it would require for comparison a more complex 3D 
stress analysis. Work in this direction is in progress. 
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A Goodman's model 



The calculation scheme involves first solving the exact model decribed in sec- 
tion |U but with fid = 0, to obtain the corresponding interfacial (symmetric) 
pressure field po{x). The second step is to solve the same constitutive equations 
for the following new boundary conditions: 

n,.(x,0) = u,(x,0) =0 (A.l) 

a.^,{x,l)+l^dCTzz{x,l)=0 (A.2) 

(Tzz{x, 1) = for |x| > a (A. 3) 

a^zix, 1) = -po{x) for \x\ < a (A.4) 

where z = and z = 1 correspond to the locations of the base and the surface 
of the elastic film respectively and po{x) is the pressure field that results from 
the first step. Eq. flA.ip accounts for the perfect adhesion of the film to its 
rigid base, Eq. (fT2l) corresponds to Amontons' law of friction with a dynamic 
friction coefficient fid and Eq. (]A.3P insures that the surface of the film is 



traction-free outside the contact zone of width 2a. Eq. (]A.4p corresponds to 
Goodman's assumption which implies that the interfacial normal stress field 
is not affected by frictional stress, and so Po{x) from the previous step is used. 

The Pourier transform of Eqs. (lA.3p - (lA.4p yields 

/ a^f) (k) cos kx dk = for |x| > a , (A. 5) 

^0 

/ afhk) coskx dk = —pix) for Ixl < a . (A. 6) 

Eq. (1A.5P is identically satisfied by 

a^hk) = r F{t)Mkt)dt (A.7) 
Jo 



where Jq{x) is the Bessel function of the first kind. By replacing (]A.7P into 



Eq. (1A.6P we get the following integral equations that determine the function 
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F(f) 

where is the function defined in Eq. ffTSjl (see section H]) obtained for the 
particular case where Hd = 0. The solution for F{x) is readily given by F{x) = 

(j){x). 

The normal stress cr^zix, 0) and the tangential stress cr^zix, 0) at the rigid base 
are then given by 

a 

(tUx,0) = J [Zi{x,t)+fidZ3{x,t)](f){t)dt , (A.9) 



a 

a^z{x, 0) = - 1 [l^dZr,{x, t) + ^3(0;, t)] (l){t)dt , (A.IO) 


where the kernels Zi{x, t) and A{k), B{k and C{k) are given by Eqs. (!27j) - (!35|) . 
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